Spherically symmetric solution in a space-time with torsion 
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By using the method of group analysis, we obtain a new exact evolving and spherically symmetric 
solution of the Einstein-Cartan equations of motion, corresponding to a space-time threaded with 
a three-form Kalb-Ramond field strength. The solution describes in its more generic form, a space- 
time which scalar curvature vanishes for large distances and for large time. In static conditions, it 
■ reduces to a classical wormhole solution and to a exact solution with a localized scalar field and 

a torsion kink, already reported in literature. In the process we have found evidence towards the 
construction of more new solutions. 
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I. INTRODUCTION 



Unification of all fundamental forces in nature has represented a difficult task along the last few decades. The most 
promising models are encoded within the context of string theory. For instance, there are models in heterotic string 
theories in which grand unifications theories are present, while in type II superstring theories, partial unification is an 
, available option in several mo del- dependent scenarios, where there exists the possibility to study geometric aspects 
' of general relativity at the microscopic level (see [l[ and references therein) . 

In order to construct effective four-dimensional theories with a minimal (non-zero) number of supersymmetric 
t— I \ generators, critical type II superstrings are typically compactified on Calabi-Yau (CY) threefolds, which break 
7— I . three-quarters of the total supersymmetry. Roughly speaking, this follows from the fact that, in the absence of 
other supergravity fields, the ten-dimensional gravitino satisfies the constraint S 1 ^^ — Va/c" 4 = 0, implying the 
presence of two four-dimensional covariantly constant spinors. Under these fluxless conditions, it is possible to 
construct a maximally symmetric four-dimensional theory (Minkowski, de Sitter or anti-de Sitter). However, fluxless 
^ , compactification scenarios suffer from the moduli stabilization problem. 

Moduli are stabilized, once one includes fluxes in the compactification process. Nevertheless, fluxes back-react and 
their presence modifies the internal geometry, forcing us to depart from the nice and smooth CY geometry. Even 
more, the presence of fluxes contributes to the ten-dimensional gravitino variation by 

S^m = V M e A + n M e A , (1) 

where km encodes the flux contributions: it vanishes in the fluxless case. Notice however, that to preserve two 
supersymmetries in the effective four-dimensional theory, the variation of the gravitino must vanish, implying that 
the four-dimensional component of the supersymmetric parameter is not covariantly constant with respect to the 
Levi-Civita equation, but with respect to a connection with torsion. It is then clear that, the effective four-dimensional 
supergravity theories obtained from the compactification process, must be different from the theories constructed in 
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the fluxless case. Such effective theories receive the name of gauged supergravities. For more details, see 1 and 0. 

One interesting question, concerns the existence of spherically symmetric solutions in the context of gauged 
supergravities and in particular, the specific case in which a Kalb-Ramond (KR) three- form flux is extended 

in the effective four-dimensional space-time. The presence of KR fluxes implies the study of gravitational theories 
in curved backgrounds with torsion Q. Torsion appears as an antisymmetric tensorial piece in the Christoffcl 
connection symbol in Einstcin-Cartan (EC) theory. The simplest level of torsion theory can provide a classical 
background for quantum matter fields and a relationship with spin [3]-@- This is important since, an alternative 
way to study geometric aspects of general relativity at the microscopic level, is given by the EC theory [l(J, which 
is characterized by the presence of a spin angular momentum in addition to the mass. The main task to include 
torsion in space-time is to unify gravity with electromagnetism at a microscopic level, where matter formations are 
performed by elementary particles. 

Torsion can be incorporated as well into the classical solutions of Einstein equations, through a non-symmetric 
part of the KR metric, which in turns reduces to the effective antisymmetric two-index potential -B M „. The field 
strength of has in general, a gauge invariant coupling with the electromagnetic field, leading to significant effects 



on many cosmological and astrophysical phenomena 



A complete understanding of torsion phenomena in the EC theory require appropriate knowledge of the various 
solutions of the gravitational field equations in space-times under various circumstances. Several cases have been 
studied in the last few years. For instance, in [14l ] SenGrupta and Sur obtained the first known solution of 
Einstein equations in a space-time with torsion (effectively induced by the presence of a four-dimensional KR field), 
corresponding to a static and spherically symmetric solution. Later on, in [13| . an approximate asymptotic solution 
was studied, to analyze the gravitational lensing and the perihelion precession, and more over, it was shown that this 
solution represents a wormhole model for a real KR field (see also }15j|). 

There are important progress in literature about the physical implications on the possible existence of torsion 1 . 
For instance, in (l7j it was reported tight bounds on many components of the torsion tensor (mixed-symmetry, 
trace and axial components) using the high-precision data from masers and torsion pendulums. See also |18| for an 
improvement on the constraints on the axial components. 

Motivated by this, our goal in this paper is to explore the possible existence of an evolving (time-dependent) 
spherically symmetric solution of the gravitational field equations (KR space-time) in presence of torsion 2 . This means 
we shall start up with an effective four-dimensional action with a torsion component in the sense of Einstein-Cartan 
theory The torsion component is given by the presence of a three-form flux which in turn is motivated from the 
string theory perspective. In particular, following [6| we consider an effective action which comes from considering 
a type IIB string compactification with Neveu-Schwarz (NS-NS) and Ramond-Ramond (RR) three-form fluxes with 
coordinates only on the extended space-time. Effectively those fluxes are interpreted as three-form fluxes which can 
be real or complex, according to the specific flux configuration selected in the string theory action. 

Once we have specified the effective action, we use gro up analysis as a method to find an analytical solution of 
non- linearized field Einstein-KR equation systems [l9j-[2l|. The group analysis lets us know if a differential system 
admits solutions with specific symmetries such as, for example, spherical symmetry. By this method, we obtain a 
set of solutions which are invariant with respect to an admitted group of point transformations. Some important 



works about symmetry analysis, can be seen in 2l|-[23j, and some applications to Einstein equations are found in 
[2lj and [24j . Therefore, by means of this method, we find an analytical evolving and spherically symmetric solution 
of Einstein equations with torsion. 



1 We would like to thank Neil Russel for pointing out these references. 

2 Although not explored in this work, such solutions could be important in order to realize the effect of torsion on electromagnetic waves 
coming from distant galactic sources |l6l . 
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Our paper is organized as follows: In section II we focus on the action of a purely gravitational field theory in 
a background with torsion and we compute the corresponding equations of motion. They consist on two set of 
equations which must be reduced; this is done in section III, where we obtain explicitly the final system equations for 
the Einstein-KR model. Using group analysis (analysis of symmetries), in section IV we find an evolving spherically 
symmetric solution. In section V, we take some appropriate limits on our solutions and we compare them with previous 
static solutions reported in literature. Finally we show our conclusions and briefly comment on some possible physical 
implications. 



II. FIELD EQUATIONS FOR GAUGE-INVARIANT EINSTEIN-CARTAN-KALB-RAMOND COUPLING 

The presence of fluxes in string compactification yields the construction of a variety of four-dimensional effective 
actions depending on the way the fluxes are turned on. Generically, for the fluxes living on the internal space is 
possible to stabilize the moduli related to the compactification process itself. 

For fluxes with coordinates on the extended four-dimensional space-time after compactification, effective actions 
can describe a curved space-time with torsion as shown in [8f , where it was considered the presence of a Neveu-Schwarz 
- Neveu-Schwarz (NS-NS) three-form flux in a ten-dimensional type IIB superstring with a bosonic action given by 

S = J_ ( d w 
2k 10 J 

in the Einstein- frame 3 , and where r = Co + ie 
Ramond (RR) zero- form Cq. 

After compactification 4 and by considering a space-time threaded with a three-form flux H^ vpi this field induces 
a torsion component on the connection 5 , yielding an effective (3 + 1)— dimensional and gauge- invariant action, with 
an Einstein-Cartan-Kalb- Ramond (ECKR) coupling in presence of external matter fields (we shall closely follow [6]). 
The action is given by 

S =IJ ' V=fl (£ (0, T) - ^H, vX H^ x + kT^ X H^ x ) d 4 x , (3) 

where R (g, T) is the scalar curvature of the EC space-time and 7" /il/A is the connection that contains the torsion ten- 
sor. Notice that is the anti-symmetrization of the affine connection in EC space-time. As usual, H^ vri — d\nB vv \ 
is the strength of the KR field, B^ v and k = 16irG is the corresponding coupling constant. 

As shown in Q, the curvature scalar R(g,T) of the EC space-time can be written in terms of the curvature of a 
purely Riemmannian space-time R{g) and the corresponding torsion tensor T^x as 

R{g,T)=R{g) + n% v xT^ x , (4) 
where T^vX satisfies the constraint equation 

T pu x — —-Hpux. (5) 



is the complex scalar formed by the dilaton <f> and the Ramond- 



3 Where we have used that the Einstein-frame and string-frame metrics are related by G^ N = e~^^ 2 Gf IN . 

4 We follow the standard notation on which capital indices M, N,P, .. are used to denote ten-dimensional coordinates, while greek indices 
fi, v, A, ... are used to describe four-dimensional extended coordinates. 

5 As briefly commented in the introduction, the torsion component can be elucidated from a stringy point of view. By requiring the 
effective four-dimensional theory to have 4 supersymmetric generators, the presence of three-form fluxes contributes with an extra term 
in the ten-dimensional gravitino variation, establishing the presence of an effective torsion component on the four-dimensional connection. 
Different flux configurations lead to a variety of four-dimensional gauged supergravities. By considering a simple configuration given by 
a single 3-form flux extended in four-dimensions, EC action should describe the bosonic part of a specific gauged supergravity theory 
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Therefore, it is followed that the KR field strength H 3 plays the role of the spin angular momentum density, 
considered to be the source of torsion Q . 

In our case, we start from the ten-dimensional action in the presence of both, NS-NS and RR three-form fluxes 
Hmnp and Fmnp respectively. The corresponding action is given by 

where Gmnp = t Fmnp ~ Hmnp is a complex three-form. Notice that by turning off the RR field F 3 , the remnant 
is a real field. Now, after compactification and by considering the complex three-form flux to have components on the 
extended four-dimensional space-time only, a torsion term in the effective four-dimensional theory is induced. This 
is analogous to the term in the action given by Eq. ([3]), with H^„\ being substituted by G^ v \- The effective action 
reads 

S = ~[v=g(R (.9, T) - -G^xG^ + kT^G^x) d A x , (7) 
which we shall also refer to as ECKR action. 

Therefore, after making use of the constraint on the torsion element, the corresponding equations of motion are 

G>„ = Rfj, u ~^QiivR — kT^ iv , (8) 
D^ uX = -Ld^ (V^gG^) = 0, (9) 

V K 

where T^, a symmetric 2-tensor, plays the role of the energy-momentum tensor 6 , which reads 



T,a, = j Ug vp G a ^G a&p - \g^G a MG aPl ^ • (!0) 

Following the notation used in [lij ]. let us denote by hi, fi2, h 3 and /14 the covariant four independent complex 
components, G012, G013, G023 and ^123, of the complex flux C/ 3 , and by h , h 2 , h 3 and h A the corresponding complex- 
conjugated contravariant components. Notice as well that from Einstein equations, the curvature of the space-time is 
given by 

i? = -^ 3 | 2 , (11) 

where \Gz\ 2 = G^\G^ X - Therefore we expect that a solution of the Einstein equations in the presence of torsion, 
describes a space-time with a negative curvature proportional to the flux, which in turn, is the source of torsion. 

Recall that our goal is to find time-dependent exact solutions for the system of Eqs. ©. For that, we shall consider 
the most general metric with an evolving spherical symmetry, 

ds 2 = n 2 (t) [e v dt 2 - e x dr 2 - r 2 {d6 2 + sin 2 8d<p 2 )] , (12) 

where v = v(t,r), A = X(t, r) and Q — Q(t) is the conformal factor, which is finite and positive-definite over the 
domain of t. Therefore, from Eqs. ((8|) and (fT0| we obtain the following set of five differential equations, corresponding 



6 Notice that this tensor is directly related to the spin-density tensor in the field equations for the ECKR space-time 
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to the non-zero components of Einstein equations: 



Rft 2 (h 1 ^ + h 2 h 2 + h 3 h 3 ~ h 4 h 4 ) 

Rn 2 (h 1 ^ + h 2 h 2 - h 3 h 3 + h A hi) 

2RQ 2 (h 1 ^ - h 2 h 2 + h 3 h 3 + h 4 h 4 ) 



X' 1 \ 1 



v' 1 \ 1 



A n + 3 U 



- - 2- + 



z/ 2 _ jAY z/ - A' 
" ~2 2~~ + r 



+2e~ 



A + 



A 2 z>A 



(13) 

(14) 
(15) 



2Ril 2 (-h 1 ^ + h 2 h 2 + h 3 h 3 + h A h A 



+2e -w 



v >2 


v'X' 


v' — X' 
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■ \ 2 






n\ 












n } 


~ 2 n' 



A + 



A 2 vX 



(16) 



2nn 2 h 3 h 4 



x ,n 



(17) 



where an over-dot denotes a derivative with respect to the time coordinate t, a prime denotes a derivative with respect 
to r and R =|k. On the other hand, the vanishing components of Einstein equations, provide with a set of constraints 
on the h fluxes, given by 



h 4 h 2 = h A h x = h 3 h 2 = h 2 h x = h x h 3 = 0. 



(18) 



Besides the above differential equation, there is another set of equations corresponding to the dynamics of the field 
H 3 given by Eq. (J9)). These six independent equations read 



= h\ + h 1 cot 9 + h 2 3 , 

o = h\ ^' + x ' 2 





= h\ 
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+ A' 


2 N 
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+ A 


-4 
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+ A 


-4 




2 
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+ A 
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-42 



h 1 - h 3 



3' 



h 2 + h\ + h 3 cot 0, 



,3- 



i/' + A' 2 



(19) 
(20) 

(21) 
(22) 

(23) 

(24) 



where the subindex i runs over (0, 1,2,3) and denotes a derivative with respect to the corresponding (t, r, 0, 0) variables. 

The system of differential equations is reduced by Eqs. ([TS jl - p^]) . which implies a vanishing values for the fluxes 
hi and h 2 . Notice that these flux values satisfy the constraint imposed by the zero components of Einstein equations 
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given by Eq. (I18[) . More reductions can be performed by considering Eq. (|20[) and Eq. 
symmetry on the system, such that 



h% = h% = 0. 



We select an axial 



(25) 



It is important to stress out, that this condition allows the parameter A and u> to be time-dependent, leading us to 
the construction of a corresponding time-dependent solution. 



The last conditions, reduces the system of Eqs. (|19l) - ([24|) to three differential equations, given by 

= h 3 2 + h 3 cot 9, 

V + A' 2" 



i> + A il \ , o . 



-h 2 — h cot i 



The corresponding solutions for ha and /14 read 

h 3 = 
h 



ha _ u 3 (t,r) 
Q 6 r 4 sin 2 Be v ~ sin6» 
—hi _ Ui(t,r) 



n 6 r 4 sin 2 9e x 



sinf 



(26) 
(27) 
(28) 

(29) 
(30) 



where we have expressed them in terms of the functions 113 and 1*4, both depending on the parameters t and r. 

Now, by substituting Eq. and Eq. ^ into Eq. (gZJ and into Eqs. ([T5 ]) -([T7 ]l we obtain the following system 
of differential equations, 



u 3 + u 4 



-Re 8w r 4 (u\e v + M 2 e A ) 



2Re 8w r 4 (u\e v - u\e x ) = 



v + \ A .\ fv' + \' 2 N 

— + 4w\u 3 -l—j- + - ] u l . 



— + — ) - — +e v [vw -2w- w 2 } , 



-2e" 



2Ke 8, " +I/ r 4 it3W4 





A' 1 
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7 + ^2 








- + ^) 
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A+ y 


— v 


'(iz-A) 


A 






+ f ?i; 


r 





1 



Xw + 3w 2 



a 1 



u'-X' 



(31) 
(32) 
(33) 

(34) 
(35) 



where we have taken fi(t) = e 1 "'*'. We are now ready to apply the group analysis method in order to obtain an 
evolving exact solution, since our system has been reduced till we have the same number of equations and variables. 
This shall be performed in the next section. 



III. KALB-RAMOND SOLUTION USING GROUP ANALYSIS 
A. Symmetry analysis and invariant solution forms 

To begin with, let us describe briefly the group analysis method 19] we shall use to find invariant solutions of 
the system of differential equations obtained in the last section. In pedestrian words, the method consists on three 
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steps 7 . First, it is necessary to find out the Lie algebra of generators of the point-transformations group. This group 
is constructed according to a particular system of differential equations. In our case, the system of Eqs. (l3l"1) - (j3"5|) 
has a specific group which we shall consider. Second, the sub-algebras, with respect to the inner automorphisms 
group, are classified. This allows to construct some forms related to the invariant solutions. Finally, substituting the 
invariant-forms, we construct a system from which we obtain the so-called factor-system, that in our case, consists 
of an ordinary differential equations system. By solving the factor-system, we construct an analytic invariant solution. 

Therefore, we start by describing the infinitesimal generators admitted by the system. Since we have a system of 
five equations, depending on two variables, in the most general form, the infinitesimal generators are given by 

X = C(x\x 2 )^- + rf(u\ . . . , u 5 ) A (36 ) 

where x 1 = t, x 2 = r, (it , u 5 ) = (113, U4, A, 1/, w), i = 1, 2 and k = 1, . . . , 5. 



At this point, it is necessary to prolongate the group operator by extending the infinitesimal generators. In case of 
a second order system, the corresponding prolongated-generators are given by 

? -*+ff^+<s4' (37) 



£f = A(rn-^A(n, (38) 



8 - D 3 (^)- P ?Me), (39) 



du a d 2 u a 

where A = £r+P?^ + P%^ + ■■■ is operator of total derivative, p? = — , pf- = and a = 1, ... ,5, 

(3 = 1,2, j = 1,2. 

By applying the prolongate group operator to every equation in our original system, we are able to solve the overde- 
termined system of differential equations with respect to coefficients C( xl , z 2 ) an d ?7 fc (u 1 , . . . , u 5 ). Schematically, this 

is given by the equation X(F) = 0, where F = denotes the initial differential equations system. This allows to 

find out the Lie algebra satisfied by the infinitesimal generators admitted by the system of Eqs. (f5T j) - (|55)) . This Lie 
algebra, denoted L 3 , has the following basis, 

d d d d 

X ^ = < ^7 + r ^T-3 U3 ^ 3u 4 — , (40) 

at or ou 3 0U4 

x > - I < 41 > 

9 d d ,,„. 

X-3 = -7T- -4m 3 - 4u 4 -— , (42) 

aw 0U3 0U4 

which satisfy the commutator relation [JTi,^] = — X2. 

Besides, in relation to the group of the inner automorphisms, there are nonequivalent sub-algebras L3, which 
according to the classification in [23], are given by the generators 

X 2 +lX z , (43) 
Ai+7^3, (44) 
X 3 , (45) 



7 For a more extensive description of this method, see [20| and |2.'jj|. 
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where 7 is a real arbitrary parameter. The sub-algebras generated by Eq. (|43l) and Eq. (|44[) provide the invariant 
solution forms. Notice that for the sub-algebra corresponding to the generator Eq. (|45)) . there is not an associated 
invariant solution because it does not contain any independent variable. 

In order to get the invariant solution form, corresponding to the sub-algebra Eq. (|4"3"]l . it is necessary to solve the 
first order differential equations, given by 

{X 2 + jX 3 ) $ = 0, (46) 
and find all invariants $. The above differential equation implies that 

dt ^ W (47) 

■i-yt A~ft 7 

Therefore, the set of invariants is given by 

$0 = r, $1 = A, $2 = V-, $3 = U3e 47t , $4 = U4e 47t , $5 = w - jt. (48) 
Now, taking <£>i = $j(<I>o), where i = 1, . . . , 5, we obtain the first form of an invariant solution 

"3 = $ 3 Me- 47t , M 4 -$4(r)e- 47 *,«; = $5(r)+7i, A = $i(r), u = $ 2 (r). (49) 

Now, since -7— = 0, it follows that $s(r) = m = const. Henceforth, the first form of the invariant solutions is given 
or 

by 

«3(r,() = $ 3 (r)e- 47i , u 4 (r, t) = $ 4 (r) e - 47t , w{t) = m + 7 i, A(r) = $i(r), i/(r) = $ 2 (r), (50) 

where m is a real parameter and <£>i (i = 1, ... ,4) are undetermined functions. It is straightforward to notice here 
that this invariant solution corresponds to a stationary system, since both A and depend only on r. In other words, 
the evolving spherically symmetric solution is obtained from the first form of invariant solution. 



By a similar procedure, the second invariant solution form corresponding to the sub-algebra (|44[) . reads 
u 3 (z,t) = <& 3 M(^r (3+47)/2 , X(z) = $!(*), u 4 (z,t) = $ 4 (z)(tr)-( 3+47 )/ 2 , 

w(t) =m + hit 7 , viz) = $2(2), z = -. (51) 

r 

We have found two forms of invariant solutions, which allow us to construct two different kind of solutions of our 
original system of differential equations. We shall find an evolving spherically symmetric solution and a dynamical 
one. These are the next and final steps. 



B. Evolving spherically symmetric solution 

Let us consider both A and v to be time-invariant. Our aim is to find an evolving non-vacuum solution in a 
background with torsion. Making use of the set of identities solutions described in Eqs. (f50|) . (|29|) and (|30j). we 
proceed to substitute them into the system of Eqs. (l3TT) - (f35|) . in order to obtain the factor-system of the ordinary 
differential equations for the undetermined functions $j(r), which is given by 

= 2# A + (52) 
\r V2 W\ ) 

(§2 <|>2\ 

= -fcir 6 -! + -± ) + r$ + Vi - 3 7 VV> 2 - 1, (53) 
= fcir 6 (M + ?T) _ >^ + ^ _ 7 2 r 2^ _ x (54) 

V W2 1plJ 4>2 



$3 $4 V 2 

h 7 — 

ll>\W% 1P2 



n ol 4*3*4^ 

= 2fcir — — ■ 1-7-7-) (56) 
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where k x = Se 8 " 1 > 0, Vi = e~ 1 > and i(j 2 = e"® 2 > 0. 

In order to solve the system of Eqs. (j52"]) - (l56l) . first we analyze their compatibility conditions. This analysis gives 
us some long terms, so we limit ourselves to explain only the process to reach the important additional information. 
Assuming $ 4 ^ (to avoid singularities into the equations), we solve Eq. ([5r>| for <J> 3 and substitute this solution 
into the system of Eqs. (f52j) -([55 l) . Under this substitution, Eq. (|53| becomes a quadratic equation with respect to 
$4. Solving for $4, we substitute the corresponding solution into Eqs. (|52|) . ([54]) and ([55]) . such that we obtain a 
system of equations for ip± and ip2- 

The obtained system reduces Eq. (|54|) . such that we get that i/j' 2 = 0, which means that ip 2 = with £ being a 
real constant. The same result can be obtained from Eq. (j52")l . If we substitute ip 2 — £ into Eq. (f55|) . we obtain the 
following differential equation for ipi, 

= \ + \r 2 ty[) 2 + Lj_ - 1 - j^V^ r# + - 1 - 3 7 Ve) Vi + (4 7 2 r 2 £ + 3) 7 9 r a £. (57) 

This last equation has the following two solutions 

i>x,a = tfr 2 + l + %, (58) 

ipi,b = ?7 2 r 2 + l + -%, (59) 

where c\ is an arbitrary constant. Using ip 2 — £, and Eq. (|58p . we can see from Eq. ([56]) that there are two kind of 
solutions: $3 = or $4 = 0. If we choose $3 = and Eq. (|58|) or Eq. (|59]l to solve the system, then necessarily 
$4 = 0. Now, if we take Eq. (|59|) and $4 = the solution must be $3 = because we need c\ ~ 0. Therefore, there 
is a unique solution in which a non-trivial torsion source can be described. This is achieved by taking $4 = and by 
selecting Eq. (|58| . 



Choosing $4 = and substituting it into Eqs. (|52|) -(f56 | . we obtain from Eq. (f55|) . an algebraic equation for $3, 
= fcir 12 $^7 2 + k 1 r 10 $l + kir s <f>lci + €c 1 r 2 + + £ 2 7 2 cir 4 , 



which solution reads, 



fci r 4 



<fr 3 =±\/f^, (60) 



where c 2 = \J —c\ and ci is allowed to be positive or negative. Important to notice that, if C\ is negative, c 2 is a real 
positive number and therefore, the corresponding flux h 3 is also real. In this case, we have considered implicitly that 
only NS-NS three-form fluxes were considered from a ten-dimensional string compactification. On the other hand, a 
positive C\ imply the presence of complex three-form fluxes, from which we are considering string compactifications 
in the presence of both, RR and NS-NS three-form fluxes. 



Finally, the complete set of solutions, using Eqs. (l50l) . (|58|) . (|60)) . $4 = and V>2 = £ are given by 



u 3 (r,t) 


- ^£i A /i p -4(™+7t) 

3r 4 V k 




u 4 (r, t) 


= 0, 


e A(r) 


= (l + ^7V) 


e u(r) 


= c\ 


fi(t) 





(61) 
(62) 



(63) 
(64) 

On the other hand (and notice that in general u 3 (r, t) is a complex function), from Eq. (|29|) we have the relation 

h 3 h 3 = | U3 (r,i)| 2 f! 6 rV, (65) 
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which reduces, once we considered the set of functions given by Eqs. (|BT1) - ([M)) . to 

/,3^ = M! e -2(m+7t). ( 66 ) 

Since we obtained that /14 = 0, from Eq. (J6TJ) we see immediately that a non-trivial solution requires that C2 and £ 
are both different from zero. Recall that torsion components in the connection are related to the presence of a flux 
C/3. In our present case, the only component which survives is ft. 3 = £023 • The squared flux, as seen from Eq. ([66l) . 
goes as r -4 at constant time, decaying to zero at infinity, faster than an ordinary electromagnetic field. Therefore, 
the corrections expected by torsion at large distances, through the presence of a 3-form flux, are smaller than those 
given by a Maxwell classical field. Notice also that, as time runs to infinity, the solution decreases if and only if 
7 > and m > 0. Therefore, we shall consider this to be the case. In such conditions, we have found an evolving 
spherically symmetric solution of the Einstein equations in presence of torsion. Finally, let us mention that, as far as 
we know, there are no reports in literature about evolving exact solutions under these conditions. 



C. Dynamic solution, second invariant solution form 



In this section, we will analyze the ECKR dynamic solutions, which represents an evolutionary system in time. We 
shall consider both, v and A, to be functions of t and r. Using the second invariant solution form given by Eq. (|51[) 
together with Eqs. (|3"TT) - (|35I) . we obtain the following factor-system of ordinary differential equations, 

= z ( p- + p- ] (z$ 4 - $3) - 2z 2 $; + 2z$ 3 + (4 7 - 3) $3 + (1 - 4 7 ) z$ 4 , (67) 



= (% ~ z 2 ) - + £i ) z ^ kl + ^ _ 3^ 27 2 - z 2 , (68) 



ipi J \ip2 

= (|V - 7) Z1& + + |±) z^h - ^ 27 2 + ^iz 2 + 2^7 - z 2 , (69) 

= ^z 4 (V4) 2 - ^z 2 W) 2 - V^'i + V 2 2 zVi' + H 2 ^2 + 2^7) z^ + 4 7 ^ 2 Vi (70) 
2^2 20i 

+ y (02 - 0iz 2 ) 0i0 2 - (tfz 2 + 2 7 2 Vi) z^ 2 + ($ 2 2 - 4> 2 ^i) 2z 4 ^ 1 fci - 2 7 V 2 2 Vi, 

= 071^7 -^ 1 -2z 4 '»'- 3 fci$3$4, (71) 

where a prime denotes a derivative with respect to z. In this section, we use the same definitions for fci,0i and tp2 
as in the last section, but here ipi, "02 also depend on z. 

It is easy to sec that if $3 = z$4, then Eq. ([67]) is fulfilled. Now from Eq. (|7T|) we have that 

2 = ^lV> 2 7 ~ M\ 2-4-y 
4 2fci 

$ 2 = (72) 

2 /c 1 

Substituting Eqs. (|72p into Eq. (|68| and Eq. ([69|) . and considering that 7 ^ — 1, we obtain two equations for ip[ and 

_ 2 ~0i z 4 2 - Z 4 0l + 7 2 ^| - 27 2 z 2 -0i0; 2 

1 Z (Z 2 l/>1 - 0' 2 )(Z 2 V'1 - 7^2) 

,/ qV>2 z2 (! - ^1 + 70i) +7(7- l>0g f7 , 
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We can see that Eq. (fTOf is satisfied identically, taking into account Eq. ([72)1 and Eq. (1751) . In other words, any 
solution of system Eqs. (|73l) gives us a solution of the factor-system. 



If we choose the special case in which 7 = — 1, hence Eq. (|68|) is equal to Eq. (|69|) . This lets us obtain one 
differential equation for the two unknown functions V'lj ifa, which reads 



V4 v> 2 _ 2 z2 0i - 3?/^ - z 2 



"01 02 ^ 1piZ 2 + 1p2 



(74) 



Making the same process to get (j?3")) . we find that equation (1701) did not vanish as parallel to the previous case, so 
the resulting equation is, 



3z 2 \^\z 2 



U'2 



+ z 2 (0 2 0'/ " ^2) + y (02 - 



-z0 2 (30 2 0i - -01-02) - 6^201 = °) 



(75) 



we can conclude for this special case that the factor system also reduces a couple of equations for two unknown 
variables, nevertheless the complexity of those equations in the general and special case give us the possibility to get 
a dynamical solution analysis for a future work. 



D. Exact evolving spherically symmetric solution for the metric with torsion 

In this section we show the explicit metric form of the ECKR space-time for some special cases. First, we substitute 
the solutions (corresponding to our first invariant solutions) given by Eqs. (|£> (|64[) into the metric shown in Eq. 
(fH?)) . As a result, we obtain the metric for the ECKR space-time, which is given by 



dg 2 = e 2(m +1 t) 



rfi 2 

T 



dr 2 



1 + £ + il 2 r 2 



-r 2 (dd 2 + sin 2 



(76) 



and the associated scalar curvature reads 



R e 2(m+ 7 t) r 4- ( 77 ) 

We see that for large distances and time running to large values, the scalar curvature reaches a zero constant value. 
Therefore, this space-time evolves into a an asymptotically flat space-time. Notice that the parameter £ does not 
play a role in the curvature values, which can be positive or negative according to the kind of fluxes, real or complex, 
giving rise to torsion. From the string theory perspective, this corresponds to turning on only NS-NS fluxes or both, 
RR and NS-NS ones. 



It is important to point out that the curvature of this solution behaves as expected from Einstein equations as 
shown in Eq. dill) , i.e., that the curvature is produced by the presence of the torsion component which in turn 
is given by the presence of the flux /13. Therefore, in the absence of the 3- form flux /13 the metric solution would 
describe a flat space-time. The curvature must be proportional to \h^\ 2 which by Eqs. fTTj) and (|66l) . is given by Eq. flTTl) . 

Notice as well that as time runs to infinity, the curvature goes to zero even around the singularity at r — 0. The cor- 
responding space-time described by this solutions is highly curved around the singularity at r = 0. However, as times 
runs, the curvature decreases and the space-time evolves into a soft universe with a null curvature for r ^ (see FigfTJ). 

Let us consider some specific values. First, notice that in the absence of flux, i.e., for c\ = 7 = m = 0, a flat 
space-time is recovered. 
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FIG. 1: Scalar curvature as a function of r and t, with £ = 7 = m = 1, and c\ = 1. 



1. For non-trivial flux, the general solution depends on four parameters (to, 7, £, ci). One interesting case for the 
branch of solutions with m = 0, involves a static configuration, with 7 = and £ = 1, and a negative c\, in 
which the metric Eq. (|76[) reduces to the wormhole solution found by Sen Gupta et.al. in [liij . where the horizon 
is found at r = c-i- Notice as well that the torsion term given by Eq. (|66|) reduces in this case to 

h 3 h 3 = A, (78) 

which represents an static flux-configuration in the presence of a real field axion (Hodge dual of the three-form 
flux). In this case, only NS-NS fluxes are taken into account in the ten-dimensional theory. 



On the other hand, this also suggests us that our exact solution Eq. (|76|) with 7^0 and to = 0, should represent 
an evolving wormhole, although a more detailed analysis is required here 8 . For the case in which to ^ and 
holding the same conditions on the rest of the parameters we obtain a static configuration with a metric given 
by 



ds 2 = 



dt 2 



dr 2 



1 + ci/r 2 



\d6 2 



sin 



(79) 



wich scalar curvature reads 



R 



-2e 



,<T 



(80) 



Notice that this is just a conformal map of the above mentioned wormhole solution, with t = t/t;. 



The presence of a three-form flux in the effective four-dimensional space-time, has been studied as well from the point of view of string 
compactifications. A relation with extremal supersymmetric black holes is shown in [29I) ■ 
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2. For m = 7 = 0, £ = 1 but ci being positive, the corresponding flux is complex. This means that from a string 
theory point of view, both fluxes RR and NS-NS are turned on. The metric Eq. (|76|) reduces in this case to the 
one studied in 30] which is given by 

ds 2 = dt 2 l —r^dr 2 - r 2 (dQ 2 + sin 2 6d 2 c/)), (81) 

1 + C\jr z 

which corresponds to an exact solution of Einstein equations with a localized scalar field (Hodge dual to a 3-form 
flux) and a torsion kink. As it was shown in [3(|, these fields reside in an open einsteinian microcosmos. For 
to 7^ we obtain also an equivalent metric through the conformal factor e 2m . For vanishing and non- vanishing 
to, there is a naked singularity at r = 0. 



IV. CONCLUSIONS AND COMMENTS 



By using the group analysis method of differential equations, we have obtained an exact evolving solution of the 
Einstein-Cartan equations of motion, corresponding to a space-time threaded with a three-form Kalb-Ramond field 
strength. The solution, related to the first invariant form constructed by this method, depends on four parameters 
namely (to, £, 7, C2), where C2 measures the number of NS-NS flux present in the configuration. 

This general dynamical solution shows some interesting features. First of all we observe that the corresponding 
scalar curvature vanishes for large time and large distances. This limit does not depend on the number of flux we 
consider, which in the present case is only related to the the flux component £023 which can be real or complex. 
The extra components of the most generic flux Q3 do not play a role in the solution. Moreover, we see that for 
large distances and time running to large values, the scalar curvature reaches a zero constant value. The parameter 
£ does not play a role in the curvature values, which can be positive or negative according to the kind of fluxes, 
real or complex, giving rise to torsion. We also notice that as time runs to infinity, the curvature goes to zero 
even around the singularity at r = 0. The corresponding space-time described by this solutions is highly curved 
around the singularity at r — 0. However, as time runs, the curvature decreases and the space-time evolves into a 
soft universe with a null curvature for r ^ 0. Therefore, this space-time evolves into a an asymptotically flat space-time. 

Special cases are commented: For instance, for the branch of solutions with to = 0, there is one interesting 
case. The solution corresponding to the values to = 7 = 0,£ = 1 and Ci < reduces to the static worm-hole 
solution found by Sen Gupta et.al. in ((HI). For generic values of parameters but to = 0, we constructed a 
solution which might be related to a dynamical evolving of a wormhole. This is studied in a forthcoming work. 
For ci > 0, the solution reduces to the spherically symmetric exact solution reported in [3(| which corresponds to 
an exact solution of Einstein equations with a localized scalar field (Hodge dual to a 3-form flux) and a torsion 
kink. Notice that c\ can be positive or negative, according to which kind of fluxes we are considering. From 
the effective four-dimensional point of view, they are related to torsion components written in terms of real 
or complex three-form fluxes accordingly. From a string theory perspective, such fluxes correspond to type IIB 
string compactifications in the presence of NS-NS fluxes (positive ci) or NS-NS and RR three- form fluxes (negative Ci). 

On the other hand, solutions characterized by the value 7 = represent static flux configurations. For the case in 
which the rest of parameters do not vanish, the solutions are asymptotically flat. 

In the process, we have found evidence that there exist more new solutions. Roughly speaking, the group analysis 
method consists on constructing invariant forms from which one could in principle, find an exact solution of the 
original set of differential equations. The first invariant form, as said, leads us to find the above described solution. 
For the second invariant form, the situation turned out to be more complicated. However, we were able to reduce 
an original set of five differential equations, into a set of only two different differential equations. The solution is 
however difficult to find and we have left it for a future work. 
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